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We show that the dynamics of phonons in a set of trapped ions interacting with lasers is de- 
scribed by a Bose-Hubbard model whose parameters can be externally adjusted. We investigate the 
possibility of observing several quantum many-body phenomena, including (quasi) Bose-Einstein 
condensation as well as a superfluid-Mott insulator quantum phase transition. 
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Systems of ultracold bosons present a rich variety of 
fascinating phenomena, like Bose-Einstein Condensation 
(BEC) |]J, or the superfluid-Mott insulator (SI) quan- 
tum phase transition |2(. At present, there exist very 
few physical systems in which these effects can be ob- 
served. Atomic gases at ultra low temperature consti- 
tute a unique system in this context, since their physical 
parameters can be adjusted using external fields, some- 
thing which has enabled the observation of BEC Q or 
the SI transitions 0, for example. In this paper we show 
that phonons in a crystal of trapped ions interacting with 
lasers provides us with another system where all these 
phenomena can be observed in a very clean way. As 
in the case of neutral atoms, the physical parameters 
describing the phonon dynamics can be adjusted using 
lasers. Furthermore, individual addressing yields novel 
possibilities for investigating novel physical situations. 

In the set up we consider here, phonons are associated 
to the motion of the trapped ions. Coulomb interac- 
tion induces the transmission of phonons from one ion to 
another, whereas anharmonicities in the trapping poten- 
tials give rise to an effective phonon-phonon interaction. 
Thus, an ion crystal is analogous to an optical lattice 
H, whereby the ions play the role of lattice sites and 
the phonons that of the atoms. An important feature of 
ion crystals is that, due to energy conservation, phonons 
cannot be created or annihilated, i.e. the phonon num- 
ber is conserved. This is in contrast with usual solid state 
systems, where phonons are subjected to processes that 
do not conserve their number thus preventing them from 
reachingj_e.g., BEC. Furthermore, the theoretical devel- 
opment H and experimental progress d, |j| 0] expe- 
rienced by the field of trapped ion quantum computation 
during the last years can be exploited in the present con- 
text to gain access to physical observables which are not 
reachable in other systems. 

Let us consider a set of N trapped ions confined by 
external electric potentials and which move around their 
equilibrium positions. The Hamiltonian describing this 
situation can be written as H = K + Vb + Vboulj where 
K describes the kinetic energy, Vq the trapping poten- 
tial, and Vcoui the Coulomb interaction between the 
ions. We will assume that: (i) the motion of the ions 
along one particular direction, say x, is decoupled from 
the motion along the other directions; (ii) the trapping 



potential along x for each ion is practically harmonic 



with frequency u, i.e. it is given by 



where 



Xi denotes the operator corresponding to the displace- 
ment of the i— th ion from its equilibrium position; (iii) 
the displacements around the equilibrium position are 
much smaller than the distances between ions; (iv) the 
Coulomb energy is small compared to the potential en- 
ergy, i.e. (3 := e 2 / (d^muj 2 ) <C 1, where do denotes the av- 
erage separation between ions. The first requirement (i) 
allows us to ignore the motional state of the ions along the 
y and z directions. Condition (ii) allows us to associate 
phonons to each of the ions in the usual way 0: if the vi- 
brational state of the i-th ion is given by the n— th (Fock) 
excitation state of the corresponding harmonic potential 
we will say that the ion has n phonons and denote the cor- 
responding state by \n). The position operator of the ion 
can be then written in terms of creation and annihilation 
operators for the phonons, i.e. Xi cx (a,- +a|). Condition 
(iii) allows us to express the Coulomb interaction between 
ions i and j as (e 2 jd\ j)xiXj, where dij is the distance 
between the ions. It is clear that this term will induce 
hopping from phonons betweeen the ions since it contains 
terms of the form a\a,j. On the other hand, condition (iv) 
imposes that the phonon number is conserved, since the 
terms of the form didj would decrease the energy by 2oj, 
something which cannot be compensated by the Coulomb 
interaction which effectively switches off this process. Fi- 
nally, the anharmonicities of the trapping potential will 
be, in lowest order, described by terms of the form xf or 
xf. Again, only energy-conserving terms will be impor- 
tant and thus only those proportional to ajetj or a] 2 a 2 
will survive. The first one will add some small correction 
to the trapping frequency, whereas the second term can 
be associated to an effective phonon-phonon interaction. 

Thus, we have shown that the dynamics of the phonons 
in an ion crystal will contain hopping terms between dif- 
ferent ions as well as on-site phonon-phonon interac- 
tions, and therefore they will be described by a Bose- 
Hubbard model (BHM). Typically, the trapping anhar- 
monicities will be very small. However, they can be 
enhanced by using off-resonant lasers. For instance, 
one may induce repulsive phonon-phonon interactions by 
placing the ions near the maximum of a standing wave, 
something which will induce an AC-Stark shift propor- 
tional to cos(kxi) 2 s» 1 — (kxi) 2 + (l/3)(kxi) 4: , where k 
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is the wave-vector of the laser. By placing them at the 
minimim, one obtains an attractive phonon-phonon in- 
teraction. 

There are several physical set-ups which realize a BHM 
as explained above. In the following we will concentrate 
in the simplest one, which consists of ions in a linear trap 
and which gives rise to a ID system. Let us emphasize, 
however, that with ions in microtraps [llj ] or in Penning 
traps |l2| one can realize higher dimensional situations. 

In a linear trap, ions are arranged in a Coulomb chain. 
Phonons moving along the chain cannot be used in the 
way we described above since for them (3 > 1 0GI1. 
However, transverse phonons corresponding to the radial 
modes fulfill /3< 1 and thus are perfectly suited for our 
purposes. The radial phonon dispersion relation has a 
large gap of the order of u>, giving rise to the phonon- 
number conservation, and has a bandwidth of the order 
of f3u). Therefore, we take x as one of the transverse 
directions and z the trap axis. The Hamiltonian that 
describes the motion in a chain with N ions is given by: 
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vboui — 



V (zi ~ z-j) 2 + (xi - Xj) 2 + (y l - yj) 2 



where uj a , a = x,y, z, are the trapping frequencies in each 
direction, and we define (3 a as the corresponding ratios 
between Coulomb and trapping energy. If u> XtV 3> u) z the 
ions form a chain along the z axis and occupy equilibrium 
positions zf. Phonons in the x direction can be described 
approximately by: 
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^x.ii are spatial dependent shifts of the trapping fre- 
quency, and effective hopping energies, respectively. Note 
that both of them are of the order of f3 x u x . The ap- 
proximations that lead to Q are: (i) In Vcouli we keep 
only second order terms in the displacements of the ions 
around the equilibrium positions. Higher order terms are 
of the form x 4 , x 2 y 2 , zx 2 . They can be neglected assum- 
ing that xo/do, zo/do <C 1, with xq, zq the size of an 
individual ion's wave-packet, xq can be estimated by the 
size of the ground state in the radial trapping frequency 
and for typical parameters (do = 5 /ito, lo x = 10 MHz, 
also used below) we have xo/d$ » 10~ 3 . In the case of 
Zq, one has to consider the collective nature of the ax- 
ial modes, because (3 Z ^> 1 if N ^> 1. If we consider 
axial modes at a finite temperature, zq is given by the 



thermal fluctuations of the position of the ion. 
limit N ^> 1, we can estimate zZ 



In the 



~ 2mu z y/0 z log 0; hu z ' 

which means that zn /rfn w 10~ 3 fc B T/Sw z , with to z = 100 
kHz, N = 100 (see filf). (ii) We consider /3 X < 1, and 
neglect the phonon number non-conserving terms in the 
couplings of the form XiXj oc (a] + a,-)(aj + dj). 

We include the effect of a standing wave in H x , such 
that a repulsive phonon-phonon interaction is induced 
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N 



|0) l (0|cos 2 (fca: l ). 



(3) 



\0)i is the internal ground state of the ions. In the fol- 
lowing we will assume that ions stay always in |0),-, and 
expand the standing-wave in the Lamb-Dicke parameter, 
rj = kx(j/hw x : 



N 1 
H m = Fj2(lW(^ + al) 2 +-v 4 (a t +al) i + 0( V 6 )). (4) 



The fourth order contribution contains a Hubbard inter- 
action, Ua\ a 2 , with U = 2Ft] 4 (note that U < if the 
ions are placed at the minimum). The other terms in 
Eq. Q are: (i) phonon conserving terms that just give 
corrections to the trapping frequency; (ii) phonon non- 
conserving terms, that rotate with frequency lo x . The 
non-conserving contributions can be adiabatically elim- 
inated if Fr] 2 /oj x <C 1. For example, in the case of the 

second order terms, Fi] 2 (a 2 + a] ), a perturbative calcu- 
lation shows that they only give harmonic corrections of 
the form ((F V 2 ) 2 /w x )(-2aJoi - 1) + O(Fr) 2 ) 3 /u 2 x . Thus 
the contributions from non-conserving term, cither give 
corrections to the trapping frequency, or can be neglected 
when compared to U. 

The final Hamiltonian takes the form of a BHM: 



N 



(5) 



where we include in H x q the corrections from the stand- 
ing wave. Note that as long as the number of phonons is 
conserved, uj x in H x q is a global chemical potential that 
does not play any role in the description of the system. 

We discuss now the properties of the solutions of the 
non-interacting Hamiltonian, H x q. A quite unexpected 
result is that the Coulomb interaction induces the con- 
finement of the radial phonons. This is due to the fact 
that the distance between ions is larger at the sides than 
at the center of the chain, and is well described by a 
quadratic dependence on the position of the ions. Thus, 
the corrections u) x ,i in Eq. @ are smaller for the ions 
placed at the center of the chain, in such a way that the 
radial phonon field is confined (Fig. The harmonic 
phonon confinement can be estimated by means of Eq. 
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FIG. 1: Phonon trapping potential (J^J for the radial modes 
of a Coulomb chain with N —50 ions, and f3 x = 10 -2 , as a 
function of the ion position along the chain. Inset: Spectrum 
of the radial collective modes, Sl x , q , that diagonalize H x q. 



© in the limit N ^> 1. In this case, the distance between 
ions at site i satisfies Il4l: 



1 



(d(i>)/d ) 



a 



3.4, 7 = 18, (6) 



where 



i — N/2. One can use Eq. Jfjjl to describe 
qualitatively the dependence of io x ,i with the position. 
We include only Coulomb interaction between nearest- 
neighbors in order to get analytical results. The spatial 
dependent part of the non-interacting boson Hamilto- 
nian is given, in this approximation, by: 



H x / {(3 x LO x 



(7) 
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(a\a i+1 + h.c). 



In the limit of many ions and low energies, the continuum 
limit in this expression describes a one dimensional sys- 
tem of bosons trapped by the frequency uj c « (8/N)f3 x L0 x . 
The lowest collective modes in the exact spectrum show a 
linear dispersion that is well described by our estimation 
for uj c (see inset of Fig. P). Thus, we get the conclusion 
that radial phonons in ions in linear Paul traps are natu- 
rally confined by an approximate harmonic potential (see 
Fig. 0). 

Our ideas lead to the following two proposals of exper- 
iments with linear Paul traps: 

(i) Super fluid- Mott insulator transition and creation of 
a superfluid phonon state by adiabatic evolution. Hamil- 
tonian JSJ describes a BHM with the peculiarity that 
hopping terms are positive, with a range that is longer 
than the usual nearest-neighbor hopping in optical lat- 
tices. However we can understand the properties of 
our system by means of the better known model with 
nearest- neighbor hopping only 0, Il6| . Let us consider 
t = (l/2)(3 x huj x , the characteristic hopping energy. If the 
total number of phonons N p h is commensurate with the 
number of ions N, then, for values U ^> t, the ground 
state of Hamiltonian JSJ is a Mott insulator, well de- 
scribed by a product of Fock states of N p h/N in each ion 
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FIG. 2: Mean phonon number {ni) = (a\ai) along Coulomb 
chains with different number of ions, N, in the state with 
Nph = N phonons in the radial lowest mode. The width of 
the wavefunction in units of N, is oc l/\/N, in accord with 
the scaling for the phonon trapping frequency, lo c oc 1/N. 



(note that phonon confinement, Hlu c is also of order i, so 
that condition U S> t ensures a uniform phonon density). 
On the other hand, the ground state for U <C t is a super- 
fluid with all the phonons in the lowest energy level. In 
Fig. |31 we present the results of an exact numerical diag- 
onalization of the complete phonon Hamiltonian (that is, 
including also the phonon number non-conserving terms) 
for the case N p h/N = 1. The transition from the super- 
fluid to the Mott insulator, with one phonon per site, 
is evident in the evolution of the phonon density as a 
function of the interaction U. 

The properties of the BHM, allows us to propose an ex- 
perimental sequence that would lead to the observation 
of the SI quantum phase transition: (1) The ion chain 
is cooled to the state with zero radial phonons by laser 
cooling. (2) Starting with a value U ^> t, the eigenstates 
of the system are well described by Fock states localized 
at each ion. The ground state of the phonon system can 
be created by means of sequencies of blue/red-side band 
transitions, in a method that has been successfully imple- 
mented with single trapped ions (see 0). (3) The value 
of U is varied adiabatically down to a given value Uf , in 
such a way that the system remains in the ground state. 
At a given critical value Uf as t, the system undergoes a 
transition to a phonon superfluid. (4) The measurement 
of the ground state can be accomplished by the coupling 
of the transverse phonons to a given internal transition. 
One could apply, for example, a red sideband pulse with 
intensity g, for a short time t. Under such conditions, 
the probability of inducing a transition to the excited in- 
ternal state is oc X)n( sm (v / "5^)) 2 -f( n ) ~ J2 n ng 2 t 2 P(n), 
where P{n) is the probability of having n phonons. This 
method would allow us to measure the mean phonon 
number. By resolving individually the photolumines- 
cence from each ion, one could observe features of the 
BEC, or SI transition in the variations of the phonon 
density along the chain. One could also apply well known 
methods to determine P(n), or even the whole quantum 
tomography of the phonon quantum states 0, ^3 . 
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FIG. 3: Mean phonon number at each ion in the ground state 
of a Coulomb chain with N = 6 ions, and N p h = 6 phonons. 
j3 x — 0.01, so that the nearest neighbor hopping terms are t ss 
5 W~ 3 uj x . Black circles represent the phonon density without 
standing wave (U = 0), and show the confinement due to 
the phonon trapping potential. Empty circles: Mott phase 
when a standing wave is applied with Frf = Q.luJ x , n 2 — 0.1, 
and U w Q.Q2ui x > t. Squares: U = 0.01u x , Triangles: U = 

5 lcr'W 



(ii) (Quasi) Bose-Einstein Condensation by evapora- 
tive laser cooling. We propose an experiment that is akin 
to the usual BEC of cold atoms in harmonic traps. First, 
we note that techniques for cooling of trapped ions, like 
laser cooling [Iolll^ | can only be used to destroy phonons. 
The existence of the trapping phonon potential in ion 
traps allows us to propose the combination of laser cool- 
ing with the idea of evaporative cooling. A possible ex- 
perimental sequence would be as follows: (1) Start with 
a Coulomb chain after usual Doppler cooling, that is, a 
chain with a given number of phonons per site, and in- 
duce a small phonon-phonon interaction U -C t, so that 
the system remains in the weak interacting regime. (2) 
Apply laser cooling at the sides of the Coulomb chain, in 
such a way that the higher energy phonons on the top 
of the confinement potential are destroyed (evaporated). 
(3) The interaction U induces collisions that thermalize 
the phonons to a lower temperature. Several cycles of 
laser cooling / thcrmalization could be applied until the 
system is cooled below the critical temperature. Detec- 
tion of the BEC could be accomplished along the same 
lines exposed above for the case of the BHM. Note that in 
the case of Coulomb chains (ID) considered here, (quasi) 
BEC is possible in finite size systems only. 

We have shown that phonons in a system of trapped 
ions can be manipulated in such a way that they un- 
dergo BEC, or a SI transition. The main ingredients of 
our proposal are: (1) The fact that phonons can have 
a large energy gap that suppresses processes that do no 
conserve the number of phonons. (2) Phonon-phonon 
interactions (anharmonicities) can be induced by placing 
the ions in a standing-wave. (3) In the particular case 
of ions in a linear trap, radial phonons would be suitable 
for this proposal, with the advantage that the Coulomb 



interaction provides us with an approximately harmonic 
phonon confinement. 

In this work we have exposed only a few applications 
of this idea, but phonons in trapped ions could be used 
to study quantum phases with a degree of controllability 
that is not possible with cold neutral atoms. Individual 
addressing would allow us to design Hubbard Hamiltoni- 
ans with local interactions that change at will from site 
to site. Different directions of the radial modes, or dif- 
ferent internal states of the ions could play the role of 
effective spins [3] for the phonons. On the other hand, 
one could also reach the regime dominated by the repul- 
sive interaction and create, thus, a Tonks-Girardeau gas 
of phonons |l9| in a Coulomb chain (this idea was im- 
plemented recently with optical lattices |2(ij). In a very 
promising approach, 2D systems of arrays of microtraps 
111, or ions in Penning traps ^3]> could be considered, 
because phonons transverse to the crystal plane satisfy 
the conditions required by our proposal. 
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